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The aim of the paper is to propose another tool for phenomenological analyses of experimental data 
from superallowed nuclear and neutron j3 decays, and from leptonic and semileptonic decays, that 
allows the finding of the most probable numerical form of the CKM matrix, as well as the determi- 
nation of decay constants, fp, and of various form factors f+(q 2 ), by using another implementation 
of unitarity constraints. In particular this approach allows the determination of semileptonic form 
factors that is illustrated on the existing data from D — > -nlv and D — > Klv decays. 

PACS numbers: 12.15.-y, 12.15 Hh, 12.15 Ff 

I. INTRODUCTION 

Within the Standard Model (SM) the flavor physics is encoded by Kobayashi-Maskawa (KM) matrix, [l| , supposed 
to be unitary, matrix that describes the quark flavor mixing through four independent parameters: three mixing 
angles, 8ij, ij = 12, 13, 23, and one CP- violating phase, S. By consequence the experimental determination of KM 
matrix entries is essential for the validation of the SM, and for detection of new physics beyond it. 

However the determination of KM entries is not an easy problem because of two different causes. The first one is 
theoretical, namely the mixing angles are not invariant quantities, their numerical values depend on the original KM 
form, 1], or on the present day form, Q, which is not rephasing invariant, see [3. These shortcomings are harmless 
if one follows Jarlskog's solution. Starting with her first papers on KM matrix, |4|, she proposed the determination 
of the quark mixing matrix in terms of directly measurable quantities, and, in the same time, invariant quantities. 
In this context an invariant quantity is one whose numerical value does not depend on the KM matrix form, or of 
its rephaising invariance. Jarlskog provided two such invariants: the KM moduli, and the celebrated J invariant, 
Ten years later, other invariants, the angles of unitarity triangles, appeared on the scene @. After few years it was 
realized that all the measurable quantities of quark mixing matrix are expressible in terms of four independent KM 
matrix moduli, see 0] and Q. In these papers it was also shown that areas of all unitarity triangles are equivalent, 
and numerically equal to half of J invariant, and in our opinion this equality also solved the J sign problem, by 
imposing J > 0. As a conclusion we can state that the numerical values for all measured invariant quantities should 
be the same irrespective of the physical processes where they are involved. 

The second difficulty comes from the experimental side. If one wants the use of KM moduli as independent 
parameters these ones are not directly measured by experimenters. In the simplest case, that of leptonic decays, they 
measure branching ratios and provide numbers for products of the form \U qq ' \fp, where U qq ' is the corresponding KM 
matrix element, and fp is the decay constant. For meson semileptonic decays the physical observable is the differential 
decay rate, dT/dq 2 , which up to known factors, is proportional to \U qq > f+(q 2 )\ 2 , where f+(q 2 ) is a complex form factor 
and q denotes the transferred momentum between initial and final mesons. In the last case the experimental teams 
usually provide numerical values for products of the form \U qq > f+(0)\, where /+(0) is the semileptonic decay form 
factor at zero-momentum transfer. It is clear that from such measurements one cannot find two unknowns, say \U qq i\ 
and fp. This can be done if and only if one can find independent constraints on KM matrix moduli; fortunately these 
ones are provided by unitarity. Thus the main aim of the paper is to show how the unitarity property of the KM 
matrix can be transformed into a powerful tool for the determination of both matrix moduli and form factors directly 
from experimental data. 

The unitarity constraints are presented in Sec. II where they are implemented in a 'x 2 — form that depends only 
on KM matrix moduli. In Sec. Ill we present the decay formulas for superallowed + — > + nuclear and neutron 
/? decays, and those for leptonic and semileptonic decays. They depend on KM matrix moduli and specific decay 
parameters such as decay constants, fp, and form factors \f(q )|, etc, that are implemented in an other x 2 -piece. In 
Sec. IV we cite the experimental papers from which we take data that are used in our fit, and in Sec. V we present 
numerical results. The paper ends by Conclusion. 
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II. UNITARITY CONSTRAINTS 



The use of \Uij\ as independent parameters raises an important problem. This means that we have to solve 
the consistency problem between moduli and unitarity property, which all amounts to obtaining the necessary and 
sufficient conditions on the set of numbers |E7y | to represent the moduli of an exact unitary matrix. After that we have 
to find a device for applying these conditions to experimental situation where data are known modulo uncertainties. 

Both these problems were recently solved, and a procedure for recovering KM matrix elements from error affected 
data was provided in Q . These unitarity constraints say that the four independent parameters Sij = sin 9ij and cos S 
should take physical values, i.e. Sij £ (0, 1) and cos 5 £ (—1, 1), when they are computed via equations set: 

T/ 2 — r 2 r 2 v 2 — Q 2 n 2 v 2 — Q 2 
v ud — c 12 c 13> v us — *12 c 13: v ub — *13 

T/2 _ „2 2 T/2 _ 2 2 

v cb — & 23 c 13> v tb ~ c 13 c 23' 

V cd = S 12 C 23 + S 13 S 23 C 12 + 2si2Sl 3 S23Cl 2 C23 COS (5, 

V cs = C 2 12 cl 3 + sj 2 S 2 13 sl 3 ~ 2s 12 S 13 S23C 12 C23 COS S, (1) 

V td = s 13 c 12 c 23 + s 12 s 23 ~ 2s 12 S 13 S 2 3Cl2C23 COS S, 

V tl = S 12 S 13 C 23 + c 12 s 23 + 2si 2 Sl 3 S23Ci2C 23 COS S 

The above relations have been obtained by using the standard KM matrix form, 2], where Vij = |t%|, and Vij 
are KM matrix entries. In paper Q it was shown that if the independent parameters are KM matrix moduli the 
reconstruction of a unitary matrix knowing its moduli is essentially unique. By consequence in the following the used 
independent parameters in all our phenomenological analyses will be Vij moduli. Although only four of them can be 
independent the experimental data "force" us to use all the possible sets of four independent moduli, as it will be 
shown in the following, and a simple combinatorial evaluation shows that there are 57 such sets. 

Relations (TTJ are rephaising invariant, i.e. they have the same form after multiplication of all KM matrix rows and 
columns by arbitrary phases. More important is that they contain all unitarity constraints. It is easily seen that all 
the six relations such as 

V* d + Vl + V 2 b = l (2) 

are a consequence of the above relations. The relations ^ are necessary conditions for unitarity fulfilment, but not 
sufficient as we show in the following. In fact relations @ show that if they are satisfied then there exist a physical 
solution for sy. Indeed from any four independent moduli entering (TTJ one can get the mixing parameters Sjj and 
cos (5, i.e. the four independent parameters entering KM unitary matrix. For example, if V us — a, V u b = b, and V c b = c 
is one set of three independent moduli, from the first five equations {T} we find all the three mixing parameters 

The other parameter, <5, can be obtained from anyone of the last four equations. If we choose the sixth equation one 
gets 

cos 5 = (1 - & 2 )(^(1 - ft 2 ) - ^) + c 2 (« 2 + b 2 (a 2 + b 2 - 1)) 
2a6cVl -a 2 - }?y/l - b 2 - c 2 

and from the remaining relations three new cos 5 formulas similar to The above relation shows that cos S is an 
other invariant in the Jarlskog sense depending of four independent moduli, and CP-violation phase can be measured 
via relations such as ((3]). 

If we make use of the last four relations ((TJ) we get only one solution for mixing parameters and cos S. Thus 
depending on the chosen four independent moduli set the number of solutions varies between one and four. Because 
there are 57 such groups one get 165 different expressions for cos S. They take the same numerical value when are 
computed via Eqs. ([T]), if and only if all the six relations similar to Eq.^ are exactly satisfied. If the moduli matrix 
generated by four independent moduli is compatible with unitarity then cosd £ (—1,1), and outside this interval 
when the corresponding matrix is not compatible. For example if we choose V us = 2257/10 4 , V u b — 359/10 5 , V c d = 
2256/10 4 , and V c b = 415/10 4 , by using the necessary relations similar to ([2]), the corresponding moduli matrix is 

/ V9490466219 2257 359 \ 

To 5 1W To 5 " i 



\U\ = 



141 3V10526471 83 

625 10' 



3 9 v i n'-i 



y/580181 V10982 V9982648619 

W 5 2500 10 5 



(5) 
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and from (|4|) one gets cos<5 ~ 0.64088, showing that the above moduli matrix, ([5]), comes from an exact unitary 
matrix. 

If we modify the previous numerical V us value by adding to it the small quantity 3 x 10~ 4 the mixing parameters 
are still physical, only s\2 is modified by a very small quantity, and respectively all the square root entries of (|5|), 
necessary for fulfilment of all the six relations similar to <J2J) . In this case one gets cos 6 —1.42427, which shows that 
the new moduli matrix, \U\, is not compatible with unitarity, even it exactly satisfies all the six relations (J2j) - 

If one computes the J invariant one finds in the above two cases 

J 2 = 6.317 x 10~ 10 , and J 2 = -1.106 x 10~ 9 (6) 

Thus the physical conditions for unitarity compatibility are cos<5 £ (—1,1), and J 2 > 0, respectively, and from a 
theoretical point of view they are equivalent. For numerical computations the use of cos 6 formulas, like (jj), seems to 
be more efficient because of their great sensitivity to small moduli variation. 

The real physical cases are those where the central value moduli matrices, directly determined from data, or from 
a fit do not exactly satisfy relations <j2j) , but only approximately; for example for a good fit the difference could be 
10 -5 — 10~ 7 , i.e. rather small from a phenomenological point of view. In these cases the different formulas for cos 6 
provide different values, physical and unphysical, even if the mixing parameters take physical values as in previous 
example. Hence physical reality obliges us to implement the unitarity constraints 

cosSi « cosSj, i ^ j, all cos6i € (—1, 1) (7) 

into a x 2 fitting device, and our choice is 

Xl= E (e V*- 1 ) + E f E ^J- 1 ) +E( COS ^ ) " COS(50) ) 2 ' -1<cos<5«<1 (8) 

j—u,c,t \i—d,s,b J j—d,s,b \i—u,c,t J i<j 

that enforces all unitarity constraints. 

III. DECAY FORMALISM 

In this section we present the decay formalism such as it is used for the description of available experimental data, 
formalism that allows us to define a second piece of the x 2 -function by taking into account as much as possible the 
physical information. 

Information on V u d come from two important sources, superallowed + — > + nuclear ft decay, and neutron j3 
decay. Superallowed + — > + /3 decay between T = 1 analog states depends uniquely on the vector part of the weak 
interaction and, according to the conserved vector current hypothesis, its experimental ft value should be related to 
the vector coup ling constant, which is a fundamental constant and by consequence has the same value for all such 
transitions, see |l0j |. (Tlj . [13 and [3. This means that the following relation should hold 

ft = 2|GvHM f | 2 = C ° nSt (9) 

where K/(hc) 6 — 2ir 3 ftln2/(m e c 2 ) 5 , Gy is the vector coupling constant for semi-leptonic weak interactions, and Mp 
is the Fermi matrix element which in this case is equal to y/2. The ft value that characterizes any j3 transition depends 
on the total transition energy Qec, the half-life, £1/2, of the parent state, and the branching ratio for the particular 
studied transition, fllj ] . The above relation is only approximately satisfied by a restricted data set, and for this set 
one defines a "corrected" value Ft = ft(l + S R )(1 + 6ns — 8c), where S R and 6ns comprise the transition-dependent 
part of the radiative correction, while 6c depends on the details of nuclear structure. In the above formula we take 
\(*v\ — 9v^ud' wu *h 9v — 1) an d write it as 

K 

^ = TT77 (10) 

where is the transition-independent part of the radiative corrections whose last estimation given in (ll| is 

= (2.361 ±0.038)% (11) 
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Similarly for neutron /3 decay we make use of formula 



Ol + 3A*) = «2£ (12) 



see 14 j, where r„ is the neutron mean life and A = gA/gv- In our approach V u d, A]£, and A are free parameters to 
be found from fit. 

In SM the purely leptonic decay of a P meson, P — ¥ li>i, proceeds via annihilation of the quark pair to a charged 
lepton and neutrino through exchange of a virtual W boson, and the branching fraction, up to radiative corrections, 
has the form 

where Gf is the Fermi constant, Mp and mi are the P meson and I lepton masses, respectively, fp is the decay 
constant, V qq > is the modulus of the corresponding KM matrix element, and rp is P lifetime. 

The next simple decays involving Vij moduli are the semileptonic decays of heavy pseudoscalar mesons, H, into 
lighter ones, P, whose physical observable is the differential decay rate, written as 

dT(H^P£v e ) _ G 2 F V q 2 q , /2 

dq~ 2 192^M£ A ( 14 ) 

where q = pn — pp is the transferred momentum, and f+(q 2 ) is the global form factor which is a combination of the 
two form factors generated by the vector part of the weak current. When the leptons are electrons, or muons whose 
masses are low compared to mass difference tuh — mp, A(g 2 ) is the usual triangle function 

\(q 2 ) = (Mfj + M 2 - q 2 ) 2 ~ AM 2 H M 2 P (15) 

For the decay B — > Dlv the experimenters make use of an other variable, namely w = (Mp + — q 2 )/ (2MbMd). 
When the r lepton is involved the above formulas are a little bit more complicated, see [l5| . 

Usually, till now, the experimenters provided numerical values for products of the form V qq i /+(0), but we expect 
that data of the form V qq i \f+{q 2 )\ will be available in the near future for all semileptonic decays. 

The second x 2 — component which takes into account the experimental data has the form 

xi = ^(^pj (16) 

where di are the theoretical functions one wants to be found from fit, di are the numerical values that describe the 
corresponding experimental data, while at is the uncertainty associated to di . For semileptonic decays di could be of 
the form di = \ f+(q 2 )\Vki ■ In the following our fitting ^-function will be 

X 2 =xl+Xi (17) 

As it is easily seen Vij moduli enter naturally in all formulas that describe leptonic and semileptonic decays being, 
in our opinion, a strong argument for their use as fit parameters. 



IV. EXPERIMENTAL DATA 

In our analysis we used superallowed + — > + nuclear /3 decays from [l(|-[ll|], and the neutron lifetime from four 
papers: llo l. 1171 . 1181 . Ilfll . for V u d determination. We also used four values for the /3-asymmetry parameter Aq, from 
papers [20|, [2l| . [22j, [23j, and one for the electron- antineutrino correlation coefficient do, [24| . for A determination. 

V us modulus is involved in kaon and pion leptonic and semileptonic decays, but also in ratio V us /V u d, as for example 
in Marciano relation, [25fj , that we write as 

V 2 f 2 B(K -> /i^( 7 ))r 7r m 7r (l - ^|) 2 

V§ l §( 1 + Cr) = =£- (18) 
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where C r is a radiative correction stemming from both n and K hadronic structures. 

The last numerical values for the product f+ K (0)V us are given by KLOE collaboration, (26|, and by FlaviaNet 
Working Group, [13], and in fit we made use of all the (little) different f+(0)V U s values corresponding to the five 
channels. We also used results from [11], [il], [1H, [ll], [l^, [11], HI], [Hi], [^], [11], [11] that give only a "mean 
value" for the above product. 

V u b is the most poo rly determined modulus although there is much experimental information coming from decays 



B — > irlv : see |4l|, [42|, [43J, |44J, [46[, 47|, [481 ]. In this papers the experimenters have been confronted with the known 
difficulty, getting two distinct parameters, V u b and \f+^(q 2 )\, from their product measured from experimental data. 
Thus they used the form factor lattice computations to obtain V ub values depending on q 2 , see [4l[, [42j], and [48j . 
For fit we found three measurements for fsVub, [Hj], (45|, [46j], and three fenomenological determinations involving 

f? n (o)v ub , 0, 0, m. 

V cd and V cs moduli enter the leptonic decays D -> /z/, [5(|, [13) [Iljand £>+ ->■ ^, [H[, [11], [H, (HI, [57j, 
respectively, as well as the semileptonic decays D — > 7r^, and £> — » ifZi/, jH^], [6(j, [Hj]. In the last three papers one 
find V cq f+(0) values, for q — d, or s, and for the first time values on 



\f + (q 2 )\v cq = x i-^ ;; 2 K * (19) 



6l| , giving the possibility of extracting form factors directl y fr om data. The above semileptonic decays allow the mea- 
surement of the ratio Vd/®'* (0) /V cs f+ K (0) , see [60j . |6l| . |63|. [63j . [6^ . [65| . and give an independent determination 
of the ratio f+ K (0)//+ ^(0), which was considered a new independent parameter. 

Finally from semileptonic decays B Dlv and B -> D*lv, ^ ^ ^ ^ ^ 
[73], one find V c b, 5(1) and J-(l) parameters. 



V. NUMERICAL RESULTS 



Data from the above cited papers were used to define x|i the second component of full \ 2 , which has a parabolic 
form in Vh. The first component Xii ®> t na t contain all unitarity constraints, has a parabolic part, and one that 
is highly nonlinear in all Vij. Thus we had to test the stability of the expected physical values against the strong 
non-linearity implied by unitarity. Eventually the chosen method was to modify all the measured central values in 
the same sense, plus and minus, respectively, proportional to their corresponding uncertainties. 

An important assumption included in our approach is that numerical Vy values must be the same irrespective of 
the physical processes used to determine them. Accordingly the other parameters, such as decay constants, fp, form 
factors, /+(0), A, etc., that parametrize each given experiment, have been considered as independent parameters to 
be obtained from fit, by applying the usual technique to obtain their mean values and uncertainties. 

The stability tests provided sets of different moduli matrices that have been used to obtain a mean value matrix 
and its corresponding error matrix. The mean and uncertainty matrices have been computed by embedding unitary 
matrices into the double stochastic matrix set, see Q. 

The central values and uncertainties of data used in fit are those published in the above cited papers, and we 
combined the statistical and systematic uncertainties in quadrature when experimentalists provided both of them. 

The numerical values obtained from fit for decay constants, f v , fx, /b, Id, Id s i semileptonic form factors /+(0), 
A^, A, as well as the ratios fx/ and, f+ K {0)/f+ v (0) are given in TABLE I. All of them are in the expected range 
although many have big uncertainties. 

Our approach allows a "fine structure analysis" of all experiments measuring one definite quantity, such as A^, or 
/?*(0). 

For example KLOE collaboration data, (26|, and FlaviaNet Working Group data, [27|, on f + (0)V us lead to 

II K (0)kloe = 950.38 ±5.56, and 

fl K {Q)Fi ama = 955.06 ±4.31, (20) 
respectively, whose central values are a little bit different, but compatible between them at lcr, and all together provide 

fl K (0)KLOE + Fla Vl a = 952.72 ± 5.30 (21) 

f+ K (0) from Table I has a precision of 1%, and the above value obtained from ten measurements has a precision of 
0.56%. 
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TABLE I: Numerical values for decay constants fp and form factors /+(0) in MeV units, and Ap] and A 



Parameter 


Central Value 


Uncertainty 


Parameter 


Central Value 


Uncertainty 


u 

fK 

!k/U 

fB 

fo 
fn s 

f+*(0) 
/f*(0) 


131.131 

154.97 

1.1818 

222.8 

207.6 

271.0 

955.34 

214.9 


1.522 

2.17 

0.0042 

25.0 

9.8 

18.0 

9.27 

13.4 


/+*(0) 

/f A "(o) 

/+ (o)//+(o)° 

HI) 

5(1) 

AK% 

A 


653.2 

751.8 

1.171 

957.5 

1,125.3 

2.373 

-1.2686 

0.002 


19.1 

10.4 

0.049 

57.7 

40.7 

0.096 

0.0057 

0.0001 






a Data on V c df+ w /V os f+ K from papers HJ-[1|], 








TABLE II: A^ central values and uncertainties provided by fit when using data from Refs. fiol 


-El 


Ref. 


[10] 
2.373(96) 


M 

2.399(112) 


[12] 
3.361(196) 


[13] 
2.294(136) 


[11] 
2.361(38)" 






a The "constant" 


A)J and its uncertainty, from [ll[ 







In this approach one could obtain information about lepton universality in leptonic decays. Because the correspond- 
ing Jm decay constant should be same for both decays M — > [iv and M — > tv, a big difference between them could 
show a possible violation. As an example we chose M = Df meson since lattice computations provided a number for 
this decay constant, fjj s = 241(3), see [78J, with a very small error. Our results are 

/ D +^ = 265.0 ±14.0 

•/W^ = 276 - 0±20 -° ( 22 ) 

which are consistent with lepton universality. The above two numbers together with that from TABLE I completely 
disagree with that provided by lattice computations, being far away from theoretical prediction at 8a, 13a, and 10er, 
respectively, where, a — 3, is lattice uncertainty. The experimental spreading is, 246.0 < f D + < 311.0, and the minimal 
and maximal values correspond to the branching ratios obtained for B(D+ — » /i + ^ M ) = (5.15±0.63±0.20±1.29) x 10~ 3 



see [58(, and to the branching ratio B(Df — > tv) = (8.0 ± 1.3 ± 0.6)%, given by Eq.(6) in paper [57j . respectively. On 
the other hand if one computes the difference between each one of the above three values and that provided by lattice 
computation, divided by the corresponding experimental error a obtained from fit one finds the same value, 1.17, 
which shows again the consistency of experimental data. In our opinion the lattice number is much underestimated. 
In contradistinction their value for Jk/I-r = 1.189(7) is not far from the fit value fx/ fn = 1.1818(42). 

Another unexpected result concerns constancy, usually assumed in all four papers ,10]-[13], assumption that is 
not confirmed by our analysis. However our result = 2.373 ±0.096 obtained from data is in good concordance 
with the value given by relation (ITT|) . see TABLE II, but our uncertainty is 2.5a higher than the theoretical one. 
In this case A)^ spreading that results from paper [lfj is, 2.2027 < A]^ < 2.472, which corresponds to 7.1cr where 
a = 0.038% is the theoretical uncertainty. The extremal nuclei are 22 Mg and 54 Co, respectively. The spreading 
obtained from Savard et al data, [l3j], corresponds to 10.8cr, and the extremal nuclei are 74 Rb and 34 C1. However the 
difference between the central results from [l(| and [l3[ is 2a, which suggests that there is still room for computation 
improvements. The mean values and the corresponding uncertainties for all data from papers [To| - 13 1 are given in 
TABLE II. 



Our approach allows the form factors determination. The paper [61| provided for the first time results on V cq \f+(q 2 ) 
q = d, s, and in TABLE III are given our determinations. Their graphic form is shown in FIG. 1 at the end of paper. 

The precision of fp and /+(0) determinations varies from 1% for f+ v (0) to 11% for Jb- More about variability of 
the above parameters could be learnt from KM moduli matrix. Our fit result for KM central moduli values is 

/ 0.974022 0.226415 0.0042512 
V c = 0.226253 0.973323 0.0381075 | (2:-!) 
v 0.0095692 0.0371307 0.999265 
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TABLE III: 


\f+(q )| form factors from D - 


■¥ nlv and D — > Klu decays 




Bin 


utt / o \ / n_i_ \ 

fl{q 2 ) (n°e + u a ) 


fl{q 2 ) (TT-e+l/e) 


f? (q 2 ) (K-e+v e ) 


f?(q 2 ) (K°e+u e ) 


1 


0.6895 ± 0.0576 


0.7072 ± 0.0355 


0.7798 ± 0.0136 


0.7808 ± 0.0208 


2 


0.7514 ± 0.0709 


0.7735 ± 0.0400 


0.8281 ± 0.0146 


0.8096 ± 0.0239 


3 


0.8442 ± 0.0841 


0.7956 ± 0.0488 


0.8435 ± 0.0156 


0.8466 ± 0.0259 


4 


0.8926 ± 0.0974 


0.9812 ± 0.0532 


0.9113 ± 0.0198 


0.8856 ± 0.0289 


5 


1.1005 ± 0.115 


1.017 ± 0.065 


0.9945 ± 0.0229 


0.9175 ± 0.0331 


6 


1 3083 + 151 


1 101 + 084 


1 007 + 027 


1 024 + 039 


7 


1.5789 ± 0.221 


1.635 ± 0.111 


1.128 ± 0.033 


1.149 ± 0.048 


8 




1.759 ± 0.235 


1.212 ± 0.042 


1.090 ± 0.062 


9 




2.024 ± 0.301 


1.303 ± 0.066 


1.233 ± 0.092 


10 






1.561 ± 0.165 


1.476 ± 0.210 



where the digit numbers are those suggested by Mathematics rounding, when working in double precision. Its 
associated uncertainty matrix is 

/ 1.1 x 1(T 6 1.9 x 10" 5 2.1 x 10~ 5 
3.6 x 10~ 5 2.7 x 1CT 4 3.1 x 1CT 4 I (24) 



\ 3.5 x 1CT 5 2.9 x 1CT 4 3.9 x 10~ 4 



The last matrix has been obtained with the help of stabilty tests. One such matrix is (|25[) that was obtained when 
all the central measured values have been modified with plus one tenth from the corresponding uncertainty. Although 
such a modification is highly improbable from an experimental point of view, it brings to light the variation direction 
for all parameters entering the fit. 

( 0.974021 0.226332 0.0074844 \ 
V+ = 0.226114 0.973083 0.044512 (25) 
v 0.0124317 0.043391 0.998891 / 

For example the above three matrices show that V u d is precisely determined with five digits, while V us , V d, and 
V cs only with three digits. V+ matrix also shows the high V u b volatility, such that future data could lead to higher 
values for it than that given by V c matrix. In fact matrices V+ and V c show that Jb € (123.5,222.8). Although our 
value for V u b = (4.25 ± 0.02) x 10 -3 is higher than that from PDG fit, 2], it is compatible with that obtained in Ref. 

48], v ub = (4.i ± o.2 st ± Q.2 syst xio '■. 

The new data from D leptonic, [53-[53, and semileptonic decays, [59|]-[6l|, combined with the new data from 
B — > Dlv and B — > D*lv (66||-[7(|, changed the KM moduli values, in particular those from the last row and column, 
see 0, p. 150. 

A big step forward will be the measurement of q 2 dependence for products of the form |/ + (q 2 )V^f,|, where q = u, c, 
for B — > nil/, B — > Dlv and B —> D*lv decays, similar to that done for D semileptonic decays. The simplest case is 
that of \f+(Q 2 )V u k\ because there are measured data in twelve bins wich can be transformed in values for the product 
\f-t-(q 2 )V u b\, [481 ] . Such measurements will allow a more precise V u b and V c b determination, and, by consequence, it 
will provide better values for all KM matrix moduli from the the first two rows, and, perhaps, the first hints for new 
physics beyond SM, if any. 

V c matrix provides numerical values for S and angles of the standard unitarity triangle, as follows 

5 = (89.96 ± 0.36)°, a = (64.59 ± 0.27)°, , . 

7 = (89.98 ± 0.06)°, /3 = (25.49 ± 0.28)° ^ ' 

The 5 value could be interpreted as a maximal violation of CP symmetry because sin<5 fa 1. A surprising change is 
the shape of the standard unitarity triangle that is now a right triangle, since 7 ~ 90°. The Jarlskog invariant is 

J = (3.567 ±0.007) x 10" 5 (27) 

Similar results are obtained from V+ matrix even if V u b is almost twice bigger than that from V c ] they are 

5 = (90.0 ± 0.2)°, a = (54.1 ± 0.1)°, , , 

7 = (90.0 ± 0.4)°, /3 = (36.0 ± 0.2)° ^ ' 
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The above results show that S and 7 angles are independent of V u b variation, while a and j3 are moderately 
dependent. Thus experimental results on the product \f + (q 2 )V u i ) \ could lead to the determination of all the angles of 
the standard unitarity triangle. 

The small angles uncertainties show that V c and V+ matrices are higly compatible with unitarity constraints and 
all 165 cos 5 formulas provide very close each other values for all of them. 

A comparison with lattice computations from shows that f£*{0) and f+ K (0), as well as their ratio are in good 
agreement with the exprimental values obtained by us, and our uncertainties are smaller. An open problem remains 
lattice computations of fs and /+"(()) parameters. 



_+ + 
D ->)tev 



0.5 1 

q 2 (GevV) 



1.5 



D -> 7i e v 



5 i 



q 2 (GevV) 



D -> K e v 



q 2 (GevV) 





1 1 1 


1 




• D°->K-e + v 

e 








1 : 

1 I : 

I 




X 


1 



q 2 (GeV V) 



FIG. 1 \f + (q 2 )\ Form Factors from D — > itlv and D — > Klv decays 



VI. CONCLUSION 



In this paper we presented a phenomcnological tool that allows determination of KM moduli, semileptonic form 
factors and decay constants directly from experimental data. It is based on a rephaising invariant implementation of 
unitarity constraints that makes use of KM matrix moduli as fit parameters. These constraints are strong enough and 
give a consistent picture of nowadays flavor physics, and until now provide no signals for new physics beyond the SM. 
However there is a discrepency between theoretical f D + lattice computation value for the — > Iv decay constant, 
and the experimental value obtained from fit, that suggests that the numerical value is highly underestimated. 

A feature of our tool is that all measurable parameters are each other strongly correlated, a little modification of 
one of them propagates to all the other parameters, property which is a consequence of unitarity constraints. 

The new data from D leptonic and semileptonic decays, and those from B — > D(D*)lv led to a significant change of 
KM moduli and of standard unitarity triangle shape. A crucial step forward would be the measurement of B —> ttIv 
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form factors in bins of 0.5 GeV 2 that will allow a better V u b determination, and a stabilization of moduli values 
entering the first two rows. 

Taking into account that our approach gives reliable results for many parameters entering flavor physics an important 
task comes to the lattice community to improve their numerical algorithms, because detection of new physics could 
be the outcome of both experimentalists and theorists. 

Acknowledgements. We acknowledge a partial support from ANCS Contract No 15EU/06.04.2009. 
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